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ABSTRACT 

Photometric stability is a key requirement for time-resolved spectroscopic observations 
of transiting extrasolar planets. In the context of the Exoplanet Characterisation Ob- 
servatory (EChO) mission design, we here present and investigate means of translating 
space-craft pointing instabilities as well as temperature fluctuation of its optical chain 
into an overall error budget of the exoplanetary spectrum to be retrieved. Given the 
instrument specifications as of date, we investigate the magnitudes of these photomet- 
ric instabilities in the context of simulated observations of the exoplanet HD 189733b 
secondary eclipse. 

Key words: space vehicles: instruments - instrumentation: spectrographs - tech- 
niques: spectroscopic - stars:planetary systems 



1 INTRODUCTION 

The last decade has seen a surge in exoplanetary discoveries, 
850 planets confirmed (Schneider et al. ]|2011 



and 



with 

over two thousand Kepler candidates (Borucki et al. |2011 



[Batalha et al. |2013|) waiting for confirmation. With the pre- 
cise measurements of their masses and radii we have gained 
a staggering wealth of information for a plethora of targets 
and planetary types. In this large and growing consensus 
of foreign worlds, some afford us the opportunity of further 
characterisation. By studying the exoplanets' atmospheres, 
we can not only infer their chemical make-up but also con- 
strain their climates, thermodynamical processes and forma- 
tion histories. The sum total of this knowledge will allow us 
to understand planetary science and our own solar system 
in the context of a much larger picture. 

Using Hubble, Spitzer as well as ground based facilities, 
the success of exoplanetary spectroscopy of transiting (e.g. 
Beaulieu et al. | 201Q, 2Qll ]|Berta et al. |2012l [C harbonneau 



et al. 2002, 2008; Demin g et al. 1 ^2005 , 2007; Grill mair et 
al. |2QQ7, 2008, Richardson et al. |2QQ7, Snellen et al. |2008| 
■ , "Stevenson et al. ||2Q10| |S wain et 

nTnTrpT 



al^P2008: "2009a"b! 'T inetti et al. ||2007a|bH20 1Qa! iTh atte et 
aL ,2010 ; Mandeh et aL |2011| [Sing et al. 12009, .201 l||Pont 



et al. 1 20081 IBurke et al. ||2007| |2010| IDesert et al. |2011| 


Redfield et al. |2008| Waldmann et al. |2012| 2013| Crouzet| 


et al. ||2012| 


|Brogi et al. |2012) as well as non-transiting 


planets (e.g. 


Janson et al. 2010| Currie et al. 2011) has 



been remarkable in recent years. 

Following these recent successes and in the frame of 
ESA's Cosmic Vision programme, the Exoplanet Characteri- 
sation Observatory (EChO) has been considered as medium- 
sized M3 mission candidate for launch in the 2022 - 2024 
timeframe (Tinetti et al. 2012b). The current 'Phase- A 



study' space-mission concept is a 1.2 metre class telescope, 
passively cooled to ^50 K and orbiting around the second 
Lagrangian Point (L2). The current baseline for the payload 
consists of four integrated spectrographs providing continu- 
ous spectral coverage from 0.5 - IG/xm at resolution ranging 
from R ^ 300 to 30. For a detailed description of the tele- 
scope and current payload design studied by our instrument 



— ^ consortium, we refer the reader to the literature ( Puig et al 



||2012| Tinetti et al. ||2012a|b 


ISwinyard et al. ||2012| |Eccle- 


ston et al. 2012| |Reess et a 


l. 112012 


Adriani et al. | 2012| 


Zapata et al. ||2012| 


Pascale et al. 


|2012| IFocardi et al. 


2012| iTessenyi et al. 


2012|. 


In order to observe the spec- 
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trum of the extrasolar planet, the EChO mission uses the 
time-resolved spectrophotometry method to observe tran- 
siting exoplanets. The spectrum of the planet is either seen 
in transmission when the planet transits in front of the star 
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along our line of sight or in emission when the planet's ther- 
mal contribution is lost during secondary eclipse when the 
planet disappears behind its host. 

By taking consecutive short observations we follow the 
transit or secondary eclipse event of an extrasolar planet. 
This approach results in an exoplanetary lightcurve for each 
individual wavelength channel. The spectrum of the planet 
is then derived by model fitting said lightcurves and record- 
ing the eclipse depths at each wavelength bin. The ampli- 
tudes of these modulations on the mean of the exoplanetary 
lightcurve depths are unsurprisingly small. Taking the case 
of the secondary eclipse emission spectrum we typically find 
the planetary contrast to be of the order of 10~^ for hot- 
Juptiers in an orbit around a typical KO star, and a much 
less favourable 10~^ contrast for so called 'temperate Super- 
Earths' orbiting M-dwarf stars (Tessenyi et al. 2012 ). A mis- 
sion such as EChO hence needs to provide a sensitivity that 
either matches or exceeds these contrasts over the duration 
of a transit/eclipse observation lasting one (e.g. GJ1214b, 
Charbonneau et al |2 009) to tens of hours (e.g. HD80606b, 
Fossey et al. ||2Q09| . The required instrument stability over 
said time spans can be achieved by providing a high pho- 
tometric stability. There are three main factors that may 
introduce photometric variability over time and may limit 
the photometric stability. These are: 

(i) Pointing stability of the telescope: 

The la pointing jitter of the satellite is currently base- 
lined to be of the order of 10 milli-arcsec from 90s to 
lOh of continuous observatior[^ Effectively this defines the 
maximum Performance Reproducibility Error (PRE) over 
ten hours, specifying the reproducibility of the experiment. 
These pointing drifts manifest themselves in the observed 
data product via two mechanisms: 1) the drifting of the spec- 
trum along the spectral axis of the detector, from here on 
referred to as 'spectral jitter'; 2) the drift of the spectrum 
along the spatial direction (or 'spatial jitter'). The effect of 
pointing jitter on the observed time series manifests as non- 
Gaussian noise correlated among all detectors in all focal 
planes of the payload and is characterised by the power- 
spectrum of the telescope pointing. The amplitude of the 
resultant photometric scatter depends on the pointing jitter 
power spectrum, the PSF of the instruments, the detector 
intra-pixel response and the amplitude of the inter- pixel 
variations. 

(ii) Thermal stability of the optical-bench and mirrors: 
Thermal emission of the instrument, baselined at 45K, is 
a source of photon noise and does hence not directly con- 
tribute to the photometric stability budge beyond an achiev- 
able signal to noise ratio (SNR) of the observation. How- 
ever fluctuations in thermal emissions constitute a source of 
correlated noise in the observations and need to be main- 
tained at amplitudes small compared to the science signal 
observed. Given a ^45K black-body peaks in the far-IR, we 
are dominated by the Wien tail of the black body distribu- 
tion, resulting in steep temperature gradients and stringent 
requirements on thermal stability in the long wavelength in- 
struments of the EChO payload. Additionally to variable 
photon noise contribution, thermal variations may impact 
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dark currents and responsivities of the detector which must 
be taken into account. 

(iii) Stellar noise and other temporal noise sources: 
Whilst beyond the control of the instrument design, stellar 
noise is an important source of temporal instability in exo- 
planetary time series measurements (i Ballerini et al. [2011 ). 
This is particularly true for M dwarf host stars as well as 
many non-main sequence stars. Correction mechanisms of 
said fluctuations must and will be an integral part of the 
EChO science study but goes beyond the instrument photo- 
metric stability discussion presented here. 

Along with the achievable SNR, the photometric stabil- 
ity of the instrument is the deciding factor for the success 
of missions or facilities aiming at transit spectroscopy. 

In this paper we study the effects of the spatial/spectral 
jitter and thermal variability on the photometric stability 
budget of the instrument and telescope. 



2 ECHO AND ECHOSIM 

Given EChO is currently in its Phase-A study phase, we 
concede that a noise budget and stability analysis, such as 
this one, can only be preliminary. Here and in Pascale et al. 
(in prep.), we present methodology used for the testing and 
optimisation of the current instrument design. The study 
presented in this paper draws on the end-to-end mission 
simulator EChOSim, which is discussed in detail in Pascale 
et al. (in prep.). 

In the simulations of this paper we assume EChO to 
have a 1.2 metre diameter primary mirror off-axis telescope. 
The light beam is simultaneously fed into five spectrographs 
via dichroic beam splitters: Visible (Vis), short wavelength 
IR (SWIR), mid- wavelength IR (MWIRl and MWIR2) and 
long- wavelength IR (LWIR) instruments covering the spec- 
tral range from 0.5 - 16/im. 

In order to provide an end-to-end simulation of the ob- 
servation, EChO Sim incorporates a full simulation of the 
science payload including the telescope, as well as the as- 
trophysical scene including Zodi emissions. EChOSim uses 
realistic mirror reflectivities and estimates the instrument 
transmission function for each channel as a function of wave- 
length. This includes transmission, optical throughput and 
spatial modulation transfer function. Using tabulated emis- 
sivity values as a function of wavelength, EChOSim also es- 
timates the thermal emission spectrum of the several opti- 
cal elements of the telescope. The transmission through the 
dichroic chain is simulated, and the incoming radiation is 
split among the 5 instruments assuming realistic transmit- 
tance and reflectivity data (Pascale et al., in prep.). Disper- 
sion and detection by the focal plane array are simulated. 
The dispersion of the light over the focal plane is modelled 
by a linear dispersion law, which is related to the sampled 
spectral resolving power R. 

The full focal plane illumination due to the backgrounds 
and the science signal's point spread function (PSF) is calcu- 
lated and convolved with realistic intra-pixel variations (see 
section|5]). Photon noise, read noise and pointing jitter noise 
are calculated on a pixel by pixel basis and time series of the 
observation are generated. These time series are then anal- 



ysed and model fitted using Mandel & Agol ( 2002 ) analyt- 



ical solutions and an adaptive Metropolis-Hastings Markov 
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Figure 1. Secondary eclipse lightcurve of a hot- Jupiter type ex- 
oplanet with echpse duration of 720min ( ^Mandel &: Agol ^2002 J . 
Noise at 10 ~^ level was added. 




Figure 2. Time series of 6 orbits of a hot-Jupiter (akin to 
HD189733b, [Torres et aini^QQSt- The deep troughs a re limb- 
darkened transits ( [Mandel L Agol ||2002| [Claret ||2000| >, smaller 
troughs are secondary eclipses and sinusoidal variations are due 
to the planetary phase curve as the planetary day-side rotates in 
and out of view. White noise of the level of 10~^ was added. 



Chain Monte Carlo (M CMC) algorithm ( [Haario et al. |2001| 
2006| [Hastings |1970| ). For a more detailed list of the instru- 
ment parameters assumed in these simulations see table [l] 
and Pascale et al. (in prep.). 



3 FREQUENCY BANDS OF INTEREST 

Figure [l] shows an example of a secondary eclipse of a typ- 
ical hot- Jupiter planet. Figure [2] shows the signal observed 
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Figure 3. Power spectra of time series shown in figure [2] for dif- 
ferent orbital periods. Blue: Period = 120 days, Green = 2.21 days 
(akin to HD189733b), Red = 0.4 days. The sensitive frequency 
range extends from 1.9xl0~^ - 1.7xl0~^ Hz. 



by EChO over the duration of 6 planetary orbits of a hot- 
Jupiter. From these figures it can easily be seen that time- 
correlated noise has the greatest impact on the retrieved sci- 
ence at temporal variation frequencies comparable to those 
of the transit /eclipse event, or multiples thereof. Figure [s] 
shows the frequency domain representation of Figure |2] given 
a variety of orbital periods. Here the desired signal is con- 
tained in discrete frequencies and their respective overtones. 
Frequency ranges beyond these shown in Figure [S] can safely 
be filtered out using pass-band filters, or normalisation us- 
ing low order polynomials in the time domain, without im- 
pairing the transit morphology. Given the range of transit 
periods observed and the goal of accurate ingress and egress 
mapping, we find the "crucial frequency band" to be from 
1.9xl0~^ to l.TxlO"*^ Hz, outside of which slow moving 
trends and high-frequency noise can effectively be filtered. 
This approach of slow moving trend removal is well tested 
for Kepler data (e.g. |Gilliland et al. [2010 ). The overall crit- 
ical frequency band for EChO is determined by the longest 
observation expected and the need to Nyquist sample the 
highest expected frequencies. We hence limit our simula- 
tions of pointing jitter or temperature fluctuations to this 
frequency range of interest. 



4 SPECTRAL JITTER 

As described in the introduction, the telescope pointing jit- 
ter has the effect of translating the spectrum on the detector 
along both, the spatial and spectral directions. Whilst the 
real translation is a combination of these orthogonal compo- 
nents, we will consider the effect of both these translations 
independently of each other. With the observing approach 
of taking consecutive spectra over a given length of time, we 
obtain a flux time series for each pixel or resolution element 
AA. Should, during the course of the observation, the stellar 
spectra be shifted along the detector the same pixel will not 
observe the same stellar flux but that corresponding to the 
shifted stellar spectrum. If we are to construct a time series 
for each resolution element of the detector, we consequently 
need to re-sample the stellar spectra to a common grid. The 
ability to re-sample the individual stellar spectra to a uni- 
form wavelength grid depends on how well we can determine 
these relative shifts. 

In this section we will simulate a series of consecutively 
observed stellar spectra and shift each one of them according 
to a pointing jitter distribution derived from the Herschel 
Space Telescope. We will flt stellar absorption lines of these 
spectra and use the derived centroids to determine the accu- 
racy with which we can determine the spectral drift. This is 
translated into a flnal post-correction flux error. It is worth 
noting that inter and intra pixel variations are less critical 
in this case as the the Nyquist sampling adopted at instru- 
ment level spreads each resolution element over at least two 
adjacent pixels. 



4.1 Synthesising pointing jitter from Herschel 
observations 

In order to simulate the pointing jitter of the EChO mission 
as realistically as possible, we synthesised the pointing pat- 
terns of the Herschel Space Telescope for which real pointing 
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Table 1. Simulation parameters for EChOSim 



Parameters 


Vis 


SWIR 


MWIRl 


MWIR2 


LWIR 


Wavelength range (/im) 


0.5 - 2.0 


2.0 - 5.0 


5.0 - 8.5 


8.5 - 11.0 


11.0 - 16.0 


Resolution (R) 


330 


490 


35 


40 


40 


Linear Dispersion (/im//im) 


15840 


6000 


689 


374 


222 


Pixel size (/im) 


90 


15 


30 


30 


30 


Detector size (pixels) 


256x10 


1111x18 


63x10 


40x10 


50x10 


Slit width (pixels) 


2 


2 


2 


2 


2 


Effective focal number (F^) 


4.0 


2.08 


1.26 


1.0 


2.0 


Effective focal length (Fe//,nm) 


5150 


2680 


1620 


1290 


2600 


Aberration parameters (Kx,Ky) 


2.5483590431 


1.4383968346 


0.9256307428 


0.6923626555 


1.356142981^ 


Dichroic emission (%) 


3.0 


3.0 


3.0 


3.0 


3.0 




jitter pdf. The Weibull distribution is given by 



0.00025^ ' ' ' 

0.0509 0.051 0.0511 0.0512 0.0513 0.0514 0.0515 

Dec (deg) 

Figure 4. Herschel pointing over a three hour continuous obser- 
vation. 



information exists (Swinyard private communication). Fig- 
ure [4] shows Herschel pointings recorded for a three hour 
time period with a sampling frequency of ^1 second. The 
pointing jitter distance for Herschel, Aner, is given by the 
Pythagorean argument relative to their means: 

Aner = VRa^ + Dec^- (1) 

Assuming the pointing jitter in Ra and Dec are nor- 
mally (Gaussian) distributed, we can state that the proba- 
bility distribution function (pdf) of Aner, P{AHer), is given 
by a Rayleigh distribution. The Rayleigh distribution is the 
Pythagorean argument of two orthogonal Gaussian distribu- 
tions. However, for a real system such as Herschel, the point- 
ing jitter is described by a Gaussian and a non- Gaussian 
component. These non- Gaussian components propagate to 
P(Aifer) in the form of an increased skew and broader wings 
than those of a pure Rayleigh distribution. We can approx- 
imate P{AHer) using the more general Weibull distribution 
of which the Rayleigh distribution is a special case. The 
Weibull distribution interpolates between a Rayleigh and an 
Exponential distribution and is hence ideally suited to de- 



where r is known as the amplitude coefficient and k as the 
distribution shape coefficient. Equation |2] reduces to an Ex- 
ponential and a Rayleigh distribution for k = 1 and k = 2 
respectively. The mean and the variance of the Weibull dis- 
tribution are given by: 

fiu^bi = tF(1 + 1/r) (3) 

^Iti = r'r(l + 2/r) - /i' (4) 
where F is the gamma function ([Riley et al. [2002 ). 



4-1.1 Scaling the Herschel pointing distribution 

Figure [5] shows P^Auer) (blue) obtained from the pointing 
information in figure [4] and the best fitting Weibull distri- 
bution, Pwbi{-), with T = 247.71 and k = 1.6547. The mean 
of the distribution is /i = 221 A5 mas. In order to obtain the 
scaled distribution to a jitter amplitude of 10 mas predicted 
for EChO, we maintain the shape parameter at At = 1.6547 
and re-derive the scaling parameter, r, for /i = 10 mas. 
This yields the pointing jitter distribution P{AEcho) shown 
in figure |6] 

We now randomly sample from P{AEcho) to obtain 
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Figure 6. Synthesised probability distribution for EChO point- 
ing jitter. 
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Figure 7. Stellar spectrum between 1.0 and 2.5 /im. Red: down 
sampled spectrum at R = 300; Blue: interpolated spectrum (see 
step 4); Inset: spectral line used to fit for the wavelength jitter; 
Gl: fiux gradient of stellar 'black-body' at 1.29 - 1.40 /xm; G2: 
fiux gradient at 1.46 - 1.64 /im. 



^Echo for M number of observations in our simulated ob- 
serving run. Assuming a resolving power of = 300 for 
the NIR channel (1.0-2.5/im), we can calculate the spec- 
tral resolving power to be 0.0159nm/mas. Hence, we can 
express AEcho as function of spectral wavelength drift and 
from hence forth denote the EChO jitter as Am (A), where 
m is the m'th spectrum observed. 



4.2 Applying pointing jitter to stellar spectra 

Using the Phoeni:5|^code we generated a stellar spectrum of 
a sun analogue and proceeded with the following steps: 

(i) The input spectrum was trimmed to a wavelength range 
of 1.0 - 2.5 /xm and re-sampled to a resolution of R = 300 at 
a central wavelength Ac = 1.75 /xm. This yields a wavelength 
coverage of ^A = 5.8 nm per pixel. We denote the spectrum 
by Fstar(A), where A is the wavelength. 

(ii) A single stellar absorption line was now selected (wave- 
length range: 1.19 - 1.21 /im), and the spectrum conse- 
quently trimmed. We denote the trimmed spectrum as 
Fiine{^) encompassing N data points in A. 

(iii) Normally distributed random noise was added to the 
spectrum with the signal-to-noise (SNR) chosen to be the 
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The resulting spectrum, Ffinai{\), is hence Ffinai{\) = 
i^Hne(A) +A/'(cr, /i), where A/'(cr, /i) is the Gaussian noise with 
cr variance and /x = mean. We chose a SNR of 500 for this 
study. 

(iv) The spectrum F/^r^a^(A) was now interpolated by a fac- 
tor S using a 'piecewise cubic Hermite interpolation' (Press 
et al. |[2007| ) which assures the best fit to the original time 
series. Here we chose S = 100. This step is required to nu- 
merically implement sub-pixel drifts and does not impair or 
bias the results. 

(v) Steps iii & iv were repeated M times to create the 
M X N dimensional matrix X, where N is the num- 
ber of points in A and M was taken to be 100. One 
can think of M being the time axis containing spectra 
Ffinai{l, X), Ffinai{2, X) . . . Ffinai{M, X), and N being the 
spectral axis of the matrix X. 

(vi) Each spectrum in X was now shifted by Am (A) along 
the wavelength axis randomly towards either the blue or red 
part of the spectrum. 



4.3 Fitting 'jittered' spectra 

(i) Now, each spectrum Ffinai{rn, A) was fitted with a Voigt 
profile using a Levenberg-Marquardt minimisation algo- 
rithm ( Press et al. |2007| and the resulting centroids. Cm (A), 
were recorded. 

(ii) The recorded centroids were subtracted from the motion 
jitter to give the fitting residual 

Rm{X) = ^(Am(A) — Cm(A))2 and the standard deviation 



of the residual an, which was recorded, 
(iii) Steps i to ii were repeated 100 times to obtain the sam- 
pled probability distribution of P{aR). 



4.4 Translating fitting residuals to total flux error 

We can now use the fitting residual Rm{X) to calculate the 
observed flux, Fobs{m, A), per wavelength range SX and time 
step m: 



Fobs (rn, A) = J F final {m, A)dA (6) 

above we can see that Rm{X) constitutes a change in the 
integration interval of each pixel with respect to the stel- 
lar spectrum. For the case where ^A 0, we can calculate 
Fobs{m, A) using the linear approximation 

Fobs (m. A) = G{Rm (A) - Rm (A)) + Ff^^ai (m. A) (7) 



G: 



dFstar{X) 

dA 



(8) 
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Figure 8. Snapshot of EChOSim output for the MWIR detec- 
tor. Top: Focal plane illumination including the science spectrum, 
instrumental and astrophysical backgrounds. Bottom: Plot of the 
total and individually contributing fluxes. 



^/2/Kx 



(10) 



(11) 



where x and y are detector coordinates along the spectral 
and spatial axes respectively, Kx and Ky are parameters 
accounting for spatial aberrations (we assume no aberrations 
Kx = Ky) and is the ratio between the effective focal 
length and the effective telescope diameter (see table [T]). 

The one dimensional PSF for the spatial direction can 
now be written as 



27T 



(12) 



Barron et al. (2007) studied the intra-pixel response of 



IR detectors, and their best-fit model is analytically imple- 
mented in EChOSim. The cross section of the response along 
the spatial axis of the detector array is hence given by: 



where ^rn(A) is the mean of the fitting residual and was sub- 
tracted to account for equal amounts of positive and nega- 
tive drifts, and G is the local gradient of the stellar spectrum. 
We can now express the resulting flux error due to residual 
drifts along the spectral direction as 



A) = G{Rm{\) - Rm{X)) 



(9) 



In section [8] we show the residual flux error Ferr for the 
gradients Gl and G2 shown in figure [71 



5 SPATIAL JITTER 

We have investigated the effect of spatial jitter as well as 
thermal variability. Whereas inter/intra pixel variations can 
largely be ignored in the case of spectral jitter, the detec- 
tor pixel responses and variations can play a significant role 
in the computation of the spatial jitter which is caused by 
the movement of the spectral PSF along the spatial direc- 
tion of the detector. EChOSim fully simulates the effect of 
inter and intra pixel variations on a spatially resolved de- 
tector grid. We here outline the spatial jitter calculations 
but refer the reader to Pascale et al. (in prep.) for a more 
detailed discussion. The dispersed signals are sampled by 
each detector assuming a wavelength-depended and instru- 
ment specific PSF which is convolved with the intra-pixel 
dependent response of the detector. Figure [S] is a snapshot 
of a standard diagnostic diagram returned by EChOSim. 
The top panel shows the focal plane illumination of the de- 
tector array. The bottom panel shows the individual flux 
contributions: astrophysical (Q point), zodiacal light (Qzodi), 
instrument thermal emissions (Q optics), and their combined 
total (Qtot)- 

We can describe the focal plane illumination of a detec- 
tor with a monochromatic point source using the diffraction 
PSF pattern (Marc Ferlet priv. com.; Pascale et al. in prep.) 



F{y) = arctan i^tanh ^ (^V + | 



arctan < tank 



2Jd 



(13) 



y - 



where 1^ is the diffusion length. The sampled PSF along the 
spatial axis is then given by the convolution 



Ps(y,X) =Py(y,X) X F(y) 



(14) 



Finally the signal is sampled by the detector. The detector 
response is given by the convolution of the point source flux 
of the star-planet system, Qpoint{X,t), with the spectrally 
dispersed PSFs 



Qpoint{iJ,t) = QE{X)Qpoint{X,t) xps{y,X) 



(15) 



where i and j are the detector pixel indices. Note that op- 
tical efficiencies don't explicitly appear in this equation as 
Qpoint{X,t) already accounts the throughput budget, with 
the exception of the quantum efficiency, QE{X) (Pascale et 
al., in prep.). 



/. ■ .\ ^ feff OPs ^ 
^P\^jJj^) — <~v ^ point 



(16) 



and cTp is the 



where feff is the plate scale (in /xm rad~"^ 
pointing jitter (in radians per second) randomly sampled 
from the pointing distribution P{AHer)- Variations in pho- 
tometric errors are estimated by consecutive runs for a range 
of pointing jitter amplitudes from zero to 200 milli-arcsec 
over a ten hour observing window. The observed spectrum 
is reconstructed using the EChOSim pipeline and the uncer- 
tainties on the reconstructed spectra are shown in figures [12] 
& [13] for two cases where PSFs of different sizes are used. 

Given the current uncertainty over the exact nature of 
inter and intra-pixel variations we have not attempted to de- 



correlate the pointing jitter in post processing (e.g. Swain 
et al. ||2008| [Burke et al~||2010| [Crouzet et al. [[2012| ) which 
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would reduce the uncertainties reported in figures p^ fc |13| 
We must hence consider this analysis as a conservative worst 
case scenario. 



6 THERMAL STABILITY 

We have also investigated the impact of thermal emissions 
and fluctuations of the optics on the photometric stability 
of the reconstructed spectrum. Whilst thermal emissions are 
not directly a problem to photometric stability (source of 
white noise) , they pose constraints on the temperature fluc- 
tuations allowed over the time span of an exoplanetary tran- 
sit/eclipse in the reddest wavelengths. Here varying thermal 
emissions may produce detector counts that are equivalent 
or exceed the signal amplitude expected from a planetary 
eclipse. We calculate the thermal contribution using 

\i + ^LAp/LD 

Qthermal{iJ,t) = J BBx^t^X (17) 

where L is the image size of the spectrometer slit in number 
of pixels, LD the lateral dispersion in mm, Ap is the pixel 
size, is the f-number and BB is the Planck function for 
a given temperature T and wavelength A. EChOSim calcu- 
lates the thermal emission of the instrument optics as well 
as the primary, secondary and tertiary mirrors separately 
and coherently propagates the resulting emission to detec- 
tor counts. We have explored a temperature regime rang- 
ing from 40 - 60K with a grid size of 0.2K for both, the 
optical bench and the mirror temperatures. This provides 
us with an absolute scale of the thermal contributions, see 
results in section [8] We have now calculated the expected 
signal strength of a secondary eclipse of a HD 189733b like 
hot- Jupiter, Fpianet, and for a given base temperature. To, 
calculated the temperature variation, AT, required to pro- 
duce a detector signal of the same strength as Fpianet for a 
given wavelength A: 



-0.002r 




No. of spectra (m) 
Figure 9. Residual centroid fitting error 




Figure 10. Probability distribution of fitting residual standard 
deviation aR 

the dynamical star-planet system and computes the time 
resolved spectra by model fitting the generated time series. 

We have now computed the emission spectra of 
HD189733b with the full noise contribution (shot noise, read 
noise, spatial/spectral jitter) as well as for the spatial jitter 
and spectral jitter only. This allows for the direct compara- 
bility of the noise contributions. 



AFpianet(X) = /a (Tq + AT) - /a(To) 



(18) 



where fx(T) is the functional form of the thermal emis- 
sion/temperature relation for a given wavelength seen in 
figures [l4] & |15| This can either be numerically approxi- 
mated by interpolation to a fine enough grid or by fitting a 
high-order polynomial function to the data. 



7 TESTCASE 

We have now calculated the spectral and spatial jitter contri- 



butions for the hot- Jupiter HD 189733b ( [Torres et al. |20Q8 ). 
We currently have insufficient information on the expected 
thermal stability in the longest wavelength range so we do 
not include thermal stability in our test-case calculations. 
This said, as seen in section [S] the overall thermal stability 
is not a decisive factor for wavelengths shorter than ~14/im. 

For our simulations we use realistic Phoenix stellar 
models for the host stars. The planetary emission model 
is taken from Tessenyi et al. (2012). EChOSim fully solves 



8 RESULTS 

As described in section |4] the 'jittered' stellar absorption line 
was fitted using a Voigt profile and the centroids recorded. 
The residual between the real and fitted spectral shifts are 



X 10"^ 



G1 X R {X) 
_G2x R^(^) 




40 60 
No. of spectra (m) 



100 



Figure 11. Fraction of residual fitting flux Ferr over the total 
flux of the star Fgtar for both spectral flux gradients shown in 
figure|7| Gl = -760 and G2 = 24 Wm-'^nm-^ . 
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Figure 12. Simulation using EChOSim showing spatial jitter 
noise as fraction of total stellar flux, studied as a function of 
spectral wavelength and pointing jitter amplitude ranging from 
200 mihi-arcseconds. A PSF FWHM is assumed of 0.7 and 0.5 of 
the detector-pixel size. Figure ^3] shows the same simulation with 
a PSF twice this size. 




Figure 15. Detector counts (in e-/s) due to the optical bench 
emission, as a function of temperature and wavelength. 
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Figure 13. Simulation using EChOSim showing spatial jitter 
noise as fraction of total stellar flux, studied as a function of 
spectral wavelength and pointing jitter amplitude ranging from 
200 mihi-arcseconds. A PSF FWHM is assumed of 1.4 and 1.0 of 
the detector-pixel size and twice the detector pixel size as stated 
in table [T] 

shown in figure [9] The Monte-Carlo analysis of the standard 
deviation of the fitting residual is shown in figure |10| As ex- 
pected, the probability distribution on the parameter gr is 
largely Gaussian as we randomly sample from the pointing 
jitter distribution P{AEcho)- This also shows that the cen- 
troid fitting does not introduce biases in the pointing jitter 
correction. The fitting residual can now be translated to a 
total flux error using equation[9] Taking the ratio Ferr/Fstar 
we can derive the relative error due to residual jitter, fig- 
ure [n] The flux error is dependent on the local stellar flux 
gradient. From figure |7| we derived two gradients: Gl = - 




Figure 14. Detector counts (in e-/s) due to the telescope emis- 
sion, as a function of temperature and wavelength. 




optics 



Figure 16. The temperature fluctuation required (ATmirror) 
for a given mirror temperature to mimic the variation of counts 
observed by a secondary eclipse of a HD 189733b like hot- Jupiter 



760 Wm- 



and G2 = +24 Wm- 



For these two 



gradients, figure [TT] shows that the relative flux error lies be- 
tween 10~^~2xl0~^. These values are of course larger for 
the Wien tail of the stellar black body wit h sp ectral jitter 



18 



& 19 



errors of ^ 10~^ at places (see figures 

Figures [12] & [13] show the simulated spatial jitter contri- 
butions for a pointing jitter amplitude range or - 200 milli- 
arcseconds. As the photometric error resulting from spatial 
jitter is largely dependent on the PSF of the instrument and 
the pixel dimensions, we have calculated the photometric 
error for a PSF FWHM of 0.7 and 0.5 of the detector-pixel 



size (figure 12) and that for double these values (figure 13). 
As seen in the figures the spatial jitter is of the order of 
^ 10~^ — 10~^ but significantly higher for the spectral ranges 
of the NIR instrument (2-5 /xm). This is to be expected as 
the SWIR instrument features a smaller pixel size. 

Figures [l4] & [15] show the thermal contribution of the 



-16 |im 
—15 |j.m 
— 14|a.m 
13|xm 
12|xm 




55 



60 



Figure 17. The temperature fluctuation required (ATopUcs) for a 
given optical bench temperature to mimic the variation of counts 
observed by a secondary eclipse of a HD 189733b like hot- Jupiter 
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Wavelength (|am) 

Figure 18. Observations of a single secondary eclipse observation 
of HD189733b binned to 0.1/im bins. Inset: native resolution of 
the instrument. 




Wavelength (|im) 

Figure 19. Error budget per wavelength in terms of planetary 
contrast ratio of HD189733b for spectrum in figure Red: total 
fiux error; blue: spectral jitter, green: spatial jitter. 

telescope (mirrors and optics) and the optical bench tem- 
peratures respectively. This can largely be regarded as neg- 
ligible for wavelengths shortward of ^14/im but significant 
at longer wavelengths and temperatures above 50K. As de- 
scribed in section |6] this thermal emission was translated 
into a temperature tolerance, ATopUcs, showing the temper- 
ature change necessary to mimic the transit depth of a sec- 
ondary eclipse feature of hot- Jupiter HD189733b at a given 
wavelength. These tolerances are shown in figures [16] &[T7| for 
the telescope and optical bench temperatures respectively. 
We find that when the telescope and optical bench are at 
45K, temperature fluctuations need to be below 6K for the 
telescope and 500mK for the bench in order to observe a 
hot Jupiter. More restrictive constraints would be required 
when fainter signals are involved. 

Having investigated the pointing jitter contributions 
for given wavelengths and pointing jitter amplitudes, we 
show the simulated 'observations' of a single eclipse event of 
HD 189733b in figure [18] The blue spectra's error-bars con- 
tain the full noise contribution. Figure [T9l shows a breakdown 
of the individual error contributions. Here red stands for the 
total noise contribution (including read and shot noise), blue 
for the spectral jitter as calculated in section ^ and green 
for the spatial jitter contribution as calculated in section |5] 



It can easily be seen that spectral jitter is important at the 
Wien tail of the stellar black body distribution and less so 
at the Ray leigh- Jeans tail. The spatial jitter noise depends 
on the individual detectors with the effect being strongest 
for the SWIR detector. 



9 DISCUSSION 

One way to reduce spatial jitter noise is to use pixels that 
are small compared to the PSF. In this case the effect of the 
spatial jitter is only governed by the intra-pixel response, 
and using many pixels to sample the PSF will washout both 
the inter and intra response variations. A higher sampling 
of the PSF using more pixels will also positively impact the 
spectral jitter as tracking of the stellar lines becomes easier. 
However, this will only be realistic if the noise from the de- 
tectors is sufficiently low but will allow a significant amount 
of de-correlation as the PSF centroid can be tracked across 
the spatial dimension of the array. The simulations take into 
account realistic intra-pixel responses in the spectral band 
from MWIR and LWIR. At shorter wavelength, the simu- 
lations currently have to be considered worst case scenarios 
because the system is not diffraction limited at these wave- 
length. 

Figures [l4] & [15] show the detected radiation in elec- 
trons/second from the instrument and telescope (separately) 
for a temperature range of 40 - 60K. Further studies will 
include additional instrumental effects such as mechanical 
vibrations, thermo-mechanical distortions, variable detec- 
tor dark currents (assumed fixed with temperature in this 
study) , detector responsitivity drifts and the effect of cosmic 
ray impacts. These studies will look at the extent the effect 
can be de-correlated from the timelines when additional in- 
formation is available for data processing, such as temper- 
ature sensors monitoring the optics and telescope tempera- 
tures. We also plan to provide off- axis detectors to monitor 
the non stellar backgrounds and therefore provide a means 
of directly removing the background signals. Such back- 
ground removal is particularly important for fainter sources. 
Whereas HD 189733b is photon limited, observations of faint 
super-earths will likely be background limited and thermal 
and zodiacal light emissions need to be carefully accounted 
for. 



10 CONCLUSION 

In this paper we present the methodology used for a pho- 
tometric stability analysis of the EChO mission and asses 
the photometric stability given its current 'Phase- A' design 
specifications. We describe how spectral and spatial jitter 
due to space-craft pointing uncertainties are propagated to 
an uncertainty on the exoplanetary spectrum measured by 
EChO. We furthermore investigate tolerances on the thermal 
stability of the space-craft's optical path. The photometric 
stability error budget was estimated for a simulated sec- 
ondary eclipse observations of the hot- Jupiter HD189733b. 
As the instrument parameters are not set in stone as of date, 
we have throughout considered the 'worst-case' assumptions 
only and photometric stability errors may significantly de- 
crease as the instrument definition phase proceeds. 
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